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Even if Supersymmetric particles are found at the Large Hadron Collider (LHC), it will be difficult 
to prove that they constitute the bulk of the Dark Matter (DM) in the Universe using LHC data 
alone. We study the complementarity of LHC and DM indirect searches, working out explicitly the 
reconstruction of the DM properties for a specific benchmark model in the coannihilation region of 
a 24-parameters supersymmetric model. Combining mock high-luminosity LHC data with present- 
day null searches for gamma-rays from dwarf galaxies with the Fermi LAT, we show that current 
Fermi LAT limits already have the capability of ruling out a spurious Wino-like solution that would 
survive using LHC data only, thus leading to the correct identification of the cosmological solution. 
We also demonstrate that upcoming Planck constraints on the reionization history will have a 
similar constraining power, and discuss the impact of a possible detection of gamma-rays from DM 
annihilation in Draco with a CTA-like experiment. Our results indicate that indirect searches can 
be strongly complementary to the LHC in identifying the DM particles, even when astrophysical 
uncertainties are taken into account. 



I. INTRODUCTION 

One of the most important questions in Dark Mat- 
ter research (DM) is to identify a possible DM candidate 
seen at the LHC with the particle responsible for the cos- 
mological relic density. Some of us have recently shown 
[l| (hereafter Paper I) that a convincing identification 
of DM particles [2|-Q can be achieved with a combined 
analysis of direct detection and accelerator data. The 
starting point of Faper I was the simulated response of 
the Large Hadron Collider (LHC) to a specific benchmark 
in the coannihilation region of a 24-parameters Minimal 
Supersymmetric Standard Model (MSSM), and the at- 
tempt to reconstruct the properties of the DM particle, 
in this case the Supersymmetric neutralino, using only 
these simulated data. As already discussed in Ref. Q, 
even 300 fb~^ of LHC data with a center of mass en- 
ergy of 14 TeV (i.e. with a dataset that should only 
become available by 2018 or so, based on current plans), 
would not allow to identify the neutralino as the sole con- 
stituent of the DM in the Universe. Assuming a standard 
expansion rate at freeze-out, in fact, the posterior prob- 
ability obtained after imposing the simulated LHC data 
exhibits multiple modes, corresponding to different neu- 
tralino compositions, and a broad peak around the true 
solution. 

We argued in Paper I that a robust and powerful way 
of breaking the degeneracy in the parameter space, and 
therefore of identifying the DM particle, is to combine 
accelerator data with ton-scale direct detection experi- 
ments measuring the recoil energy of nuclei struck by 



DM particles (see Refs. 0-111 and references therein), 
that should become available over a similar timescale. 
We demonstrated that a simple Ansatz on the local den- 
sity of DM particles (which assumes that the relative 
abundance of the particles discovered in accelerators with 
respect to DM is the same locally as their relic density 
in the Universe) is sufficient to eliminate the spurious 
mode in the posterior distribution and to constrain the 
DM properties around the true benchmark value. Simi- 
larly to the case of accelerators, a single (realistic) direct 
detection experiment can hardly provide an accurate de- 
termination of the DM properties [1, Q, but a second 
direct detection with a different target would actually al- 
low a much more precise determination of the Weakly 
Interacting Massive Particle (WIMP) mass [l3|, and if 
the new target is sensitive to the spin-dependent contri- 
bution of the WIMP-nucleus cross section it could even 
be used to discriminate among WIMP candidates [TTj . 

Here, we investigate the impact of adding information 
from indirect, instead of direct, DM searches. This de- 
tection strategy is based on the search for the annihi- 
lation or decay products of DM particles, such as high 
energy photons, neutrinos, and anti-matter 0, @- There 
are advantages and disadvantages with respect to direct 
searches. The most obvious advantage is that indirect 
searches do not require dedicated experiments. Although 
DM physics has historically played a role in establishing 
the phy sics case of experiments such as the Fermi LAT 
[13, (III, HESS [il [i| and IceCube [H, [T^, there is a 
broad range of other astrophysical motivations that made 
the construction of these experiments possible. Among 
the disadvantages, the biggest one arises from the large 



2 



uncertainties in the predicted annihilation rates (which in 
turn are a consequence of our poor knowledge of the dis- 
tribution of DM in the Galaxy and in other astrophysical 
structures), as well as in the astrophysical backgrounds 
(see also the discussion in Ref. [^). 

Gamma-rays arc often considered as ideal messengers 
for indirect detection studies, since they are not signifi- 
cantly affected by diffusion or energy losses in the local 
universe. The most promising targets are the Galactic 
center (GC) [HQ and substructures in the Milky Way 
halo, including dwarf galaxies j2J-|3l| and intermediate 
mass black holes (IMBHs) |32h39| . 

Although it has long been considered as the optimal 
target, the Galactic center is actually a very problem- 
atic region. The first big obstacle to a reliable identi- 
fication of DM, or at least for the derivation of robust 
upper limits on the annihilation flux, is represented by 
the large uncertainties on the DM distribution in a region 
which is largely dominated by baryons (iol - li^ . Further- 
more, there are strong and poorly understood astrophysi- 
cal backgrounds at the GC, that significantly complicate 
the extraction of a DM signal (see also the discussion 
in Ref. [II]). As for IMBHs, the formation scenarios 
are not very predictive and rather uncertain, therefore 
although they might be discovered as a class of objects 
with identical gamma-ray spectrum and no astrophysical 
counterparts, the non-detection does not set stringent 
constraints on the DM particles. 

Here, we focus instead on dwarf spheroidal galaxies 
(dSphs), and analyse the dramatic implications for the 
reconstruction procedure of DM properties from future 
LHC data combined with the current Fermi LAT com- 
bined analysis of 10 dwarfs [i^. We demonstrate that the 
spurious solution at low relic density is ruled out by Fermi 
LAT, if one assumes that the particle found at the LHC 
makes up all of the DM in the Universe. We also discuss 
the constraints that should become available in the next 
few years with the upcoming ground-based gamma-ray 
experiment CTA, and with a suitable analysis [47l - [5l| of 
Planck satellite data. 

The paper is organized as follows: in Section 2 we 
describe our theoretical setup and the implementation 
of experimental constraints. In Section 3 we show how 
the constraining powet of future LHC data can be im- 
proved by including current Fermi LAT constraints on 
the gamma-ray flux from dwarf galaxies, future Planck 
data probing the reionization history of the Universe and 
a possible detection of gamma-ray from Draco with CTA. 
We conclude in Section 4. 



II. THEORETICAL SETUP 

A. Benchmark Model 

For concretcness, we will start by defining our theoret- 
ical framework. In this work we will consider that the 
MSSM is the correct description for physics beyond the 



Standard Model and that it also provides a solution to 
the DM problem in terms of the lightest neutralino, Xi- 

Neutralinos are physical superpositions of the super- 
partners of the B and W gauge bosons (Bino and Wino, 
respectively) and Higgs bosons (Higgsinos) and their phe- 
nomenological properties are extremely sensitive to their 
specific composition. In particular, their annihilation 
cross section in the early Universe (which determines 
their relic abundance) as well as in DM haloes (which 
affects the gamma-ray fiux that can be observed in indi- 
rect DM detection) has an uncertainty of several orders of 
magnitude depending on the neutralino composition. It 
is worth remembering in this sense that neutralinos with 
a large Wino or Higgsino component are known to have a 
larger annihilation cross section (and hence a smaller relic 
density) than those in which the Bino component domi- 
nates. For this reason, obtaining the neutralino compo- 
sition in the LHC is a key ingredient for being able to 
determine whether or not this particle is the main com- 
ponent of the DM. 

We follow here the theoretical setup and the notation 
in Paper I: we adopt a minimal supersymmetric extension 
of the Standard Model (MSSM) with 24 free parameters, 
corresponding to its CP-conserving version. The input 
parameters are the coefficients of the trilinear terms for 
the three generations, the mass terms for gauginos (for 
which no universality assumption is made) , right-handed 
and left-handed squarks and leptons, the mass of the 
pseudoscalar Higgs, the Higgsino mass parameter /x, and 
finally the ratio between the vacuum expectation values 
of the two Higgs bosons tan /?. We assume a benchmark 
point corresponding to the low-energy extrapolation of 
model LCC3 as defined in Ref. [3]. This benchmark is 
representative of SUSY models in the coannihilation re- 
gion, where the lightest neutralino is almost degenerate 
in mass with the lightest stau. In this region, coannihila- 
tion effects reduce the neutralino relic abundance down 
to values compatible with the results from the WMAP 
satellite [S^, and therefore, the mass difference between 
the neutralino and the lightest stau is a fundamental pa- 
rameter for the reconstruction of the relic density. 

For this MSSM model, the structure of the neutralino 
sector is very characteristic of the MSSM, with the light- 
est neutralino is Bino-like, the second lightest neutralino 
is Wino-like and the two heavier neutralinos correspond 
to Higgsino-like states with a relatively large mass (ap- 
proximately 460 GcV) . This structure is a consequence of 
the low-energy hierarchy Mi < M2 < fJ. among the gaug- 
ino and Higgsino mass parameters at low-energy. Be- 
cause of this, the lightest chargino is also Wino-like, the 
heaviest corresponding to the charged Higgsino. 



B. Statistical analysis 

In order to constrain the parameters x described above 
of our 24-dimensional SUSY model, we make use of 
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Bayes' theorem 



p(x|d) 



p(d|x)p(x) 
p(d) 



(1) 



which updates the so-called prior probability density 
p(x), representing the (lack of) knowledge on the 24- 
dimensional space before taking into account the experi- 
mental constraints, d, to the posterior probability func- 
tion (pdf) p(x|d). The latter is the probability density 
after the data have been taken into account via the likeli- 
hood function, p(d|x) = £(x). Furthermore, on the RHS 
of Eq. ([T]), p{d) is the Bayesian evidence which, in our 
case, can be dropped since it simply plays the role of a 
normalization constant for the posterior in this context 
(see Ref. 53| for further details). 

The posterior encodes both the information contained 
in the priors and in the experimental constraints, but, 
ideally, it should be largely independent of the choice of 
priors, so that the posterior inference is dominated by the 
data contained in the likelihood. If some residual depen- 
dence on the prior p(x) remains this should be consid- 
ered as a sign that the experimental data employed are 
not constraining enough to override completely different 
plausible prior choices and therefore the resulting poste- 
rior should be interpreted with some care, as it might de- 
pend on the prior assumptions [s^ . [ssf . For the practical 
implementation of the Bayesian analysis sketched above 
we employed the SuperBayeS code [H, HI], extending 
the publicly available version 1.5 to the 24 dimensions 
of our SUSY parameter space. To scan in an efficient 
way the SUSY parameter space we have upgraded the 
MultiNest [HI] algorithm included in SuperBayeS to the 
latest MultiNest release (version 2.7), with uniform pri- 
ors on the pMSSM parameters as in Paper I, to which 
the reader is referred for full details. 

The set of LHC simulated measurements that we use 
as constraints in our analysis corresponds to that in Ta- 
bic 6 of Ref. Q , which assumes an integrated luminosity 
of 300 fb~^. Furthermore, as pointed out in Ref. [s^, the 
neutralino-stau mass difference can be measured with an 
accuracy of 20% with a luminosity of 10 fb~^ in models 
where the squark masses are much larger than those of 
the lightest chargino and second-lightest neutralino, as 
is our case. We therefore also include a measurement 
of the neutralino-stau mass difference in our likelihood. 
We summarize in Table U the set of LHC measurements 
adopted in our likelihood. Each of the constraints listed 
in TableUis implemented in the likelihood as an indepen- 
dent Gaussian distributed measurement around the true 
value /ii for that observable, with standard deviation ct^, 
as given in Table U i.e. the likelihood from the LHC has 
the form: 



-Clhc (x) oc cxp 



2 



M»(x))' 



(2) 



Mass 




Benchmark value, /i 


LHC error, a 


mi^i) 




139.3 


14.0 


m(x§) 




269.4 


41.0 


m(ei) 




257.3 


50.0 


m(jli) 




257.2 


50.0 


m{h) 




118.50 


0.25 


m{A) 




432.4 


1.5 


m(fi) - 


m(x?) 


16.4 


2.0 


m{uR) 




859.4 


78.0 


m{dR) 




882.5 


78.0 


m{sR) 




882.5 


78.0 


m[cR) 




859.4 


78.0 


m{uL) 




876.6 


121.0 


m{dL) 




884.6 


121.0 


misL) 




884.6 


121.0 


m(cL) 




876.6 


121.0 


m{bi) 




745.1 


35.0 


771(62) 




800.7 


74.0 


m(ti) 




624.9 


315.0 






894.6 


171.0 


771(62) 




328.9 


50.0 






328.8 


50.0 



where ^i(x) is the predicted value for the observable at 
point X in parameter space. 



TABLE I. Sparticle spectrum (in GeV) for our benchmark 
SUSY point and relative estimated measurements errors at 
the LHC (standard deviation a). 



After imposing these constraints, the ID marginalized 
posterior pdf for the neutralino relic abundance is shown 
in Fig [1] Despite the good accuracy in the determina- 
tion of some of the supersymmetric masses (for example, 
the neutralino mass has a relatively small 10% uncer- 
tainty) and especially the mass difference between the 
neutralino and the stau next-to-lightest supersymmetric 
particle (crucial to quantify the coannihilation effect), the 
resulting prediction for the relic abundance spans four or- 
ders of magnitude. This result clearly illustrates how the 
LHC might be unable to quantify the neutralino relic 
abundance and thus to determine whether or not it is 
the main ingredient of the DM. Furthermore, the poste- 
rior pdf shows a multimodal structure, with two fairly 
well separated probability density peaks indicating two 
physically different solutions. 

This uncertainty in the neutralino relic density is a di- 
rect consequence of the impossibility of determining its 
composition in an unambiguous way 0. Only the two 
lightest neutralino mass cigenstatcs arc measured (and 
none of the charginos), which is not enough to constrain 
the neutralino mass matrix. We illustrate this in Fig. [51 
where the posterior pdfs for the parameters in the neu- 
tralino mass matrix (Mi, M2, and n) obtained with LHC 
only data are shown by the empty contours. As we can 
observe, there are two possible solutions that satisfy the 
LHC constraints: one in which the neutralino is mostly 
Bino-like {Mi < M2 < /i) and another one in which it 
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depend on the neutralino self-annihilation cross section, 
{<7v), on the neutralino mass, m^o and on the spectrum 
of standard model particles produced in the annihilation 
of neutralinos, which enters in the calculation of the to- 
tal photon spectrum per annihilation (relevant for the 
search of gamma-rays from dwarfs) and of the fraction of 
energy that couples with the gas during recombination 
(relevant for CMB constraints). Indirect searches can 
therefore be used to constrain these parameters, under 
specific assumptions on astrophysical quantities, and as 
we shall see, they allow to exclude portions of the phe- 
nomenological parameter space that would remain viable 
under future high-luminosity LHC measurements. 



Fermi LAT constraints from dwarf Spheroidal 
galaxies 



FIG. 1. 1-D marginalized posterior probability density func- 
tion (pdf) for the neutralino relic abundance after LHC hypo- 
thetical measurements (given in Table 1) are taken as experi- 
mental constraints. The true value is given by the yellow/red 
diamond. 



is mostly Wino-like (A/2 < Mi < /x). Moreover, the /i 
parameter is not well determined (since the heavier mass 
eigenstates are not measured) and varies in a wide rangtQ. 
This implies that the Higgsino composition of the lightest 
neutralino can vary significantly. 

Prom the discussion above it is easy to identify which 
is the neutralino composition associated to the different 
solutions of the relic density in the pdf. The peak with 
ri-^o/i^ > 0.1 in Fig. [T] corresponds to points in the pa- 
rameter space in which the neutralino is Bino-like (thus 
having a smaller annihilation cross section and, conse- 
quently, a larger relic abundance, compatible with the 
relic density measured with cosmological data) . The long 
tail originating in this maximum and extending towards 
smaller values of fl^oh^ is obtained for neutralinos with 
an increasing value of the Higgsino component (this is, 
those for which the fj, parameter is smaller in Fig. [2]). 
Finally, the second peak situated at ^^oh"^ w 10"'^ cor- 
responds to Wino-like neutralinos, which annihilate very 
efficiently in the early Universe. 



III. INDIRECT DETECTION CONSTRAINTS 

The indirect detection signals discussed here (gamma- 
rays from dwarfs and modifications of CMB spectrum) 



^ Notice that although in theory we could have also obtained a so- 
lution in which the second-lightest neutralino was Higgsino-like 
(that is, fj, < M2), this possibility is constrained by the determi- 
nation of masses in the Higgs sector, since a light pseudoscalar 
would have also been present. 



The first 24 months of data obtained by the Fermi LAT 
telescope in survey mode have been analyzed in Ref. [4^ 
to search for gamma-ray emission from the position of 
10 dShps including Draco. The lack of detection allowed 
to set constraints on the gamma-ray emission from each 
dSph and, assuming a certain DM content, on the DM 
self-annihilation cross section (crv). The gamma-ray flux 
at energy due to DM annihilations from the direction 
^' is given by: 



1 (av) dN^ 



47r2m?o dE. 

Xi 



7 J An 



p2(r(s,*)) 



ds dfl 



(3) 

where dN^/dE^ ~ dN^ /dE^ is the total differential 
photon spectrum per annihilation, obtained by adding up 
the contributions of all annihilation channels /, weighted 
by the corresponding branching ratio Bf. The DM distri- 
bution in the dwarf galaxies p{r) is assumed to be spher- 
ically symmetric, and therefore it is a function only of 
the radius r, which can itself be expressed as a function 
of the distance along the line-of-sight from the observer 
s, and the angle with respect to the center of the dwarf 
^f. To obtain the annihilation rate, the square of the DM 
density is then integrated along the line-of-sight s over 
the solid angle A^l. 

If no excess emission is detected from the direction of 
a dwarf (which is the case so far), then Eq. ^ allows one 
to translate an upper limit in flux into an upper limit on 
the DM parameters, once a specific DM profile is assumed 
for the dSph. In Ref. [H], the Fermi LAT collaboration 
combined the data from 10 dSphs into a single likelihood 
analysis, obtaining an upper limit on (erf) of the order 
of 10^^'''cni'^s~^ for a DM mass around 130 GeV (in the 
case of annihilation into bb). The analysis in Ref. [i^ 
accounts for the astrophysical uncertainties on the DM 
profile for each dSph. The DM profile can be determined 
from kinematic data of the member stars, and in partic- 
ular measurements of stellar velocity dispersion can be 
used to build a likelihood function that depends on the 
parameters defining a DM halo profile |3l| . These quan- 
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FIG. 2. 2-Dimensional marginalized posterior pdf in the (Mi, M2) and in the (M2, /i) planes. The inner and outer contour 
encloses 68% and 98% probability regions, respectively. The empty contours are for the case where LHC only data are applied, 
whereas the filled regions include current Fermi LAT upper limits from the combined analysis of dSphs. The best fit is shown 
with the encircled black cross while the true value is given by the yellow/red diamond. 



titles have then been included in the likelihood analysis 
of Fermi LAT data, so that their final result accounts for 
our relatively poor knowledge of DM in dSphs (see Ref. 
[46| for more details). 

We focus here for definitiveness on the upper limit on 
{<jv) derived for a DM particle annihilating to b quarks, 
taken from Fig. 2 of Ref. [1^. That same figures also 
indicates how the upper limit depends on the dominant 
annihilation channel. With a neutralino mass around 
130 GeV, as is the case for our benchmark point, the 
constraints are all within a factor of two from the case 
of annihilation into b quarks, and this is not enough to 
change our results significantly (see later). The only ex- 
ception is for an annihilation predominantly into /i'^/i", 
for which the Fermi LAT upper limit is approximately 
an order of magnitude weaker. However, we checked that 
the branching ratio into muons is subdominant (with a 
branching ratio smaller than about 0.1) for all the sam- 
ples in our scan, once LHC data are included, hence this 
case can be discounted. 

We include the information provided by Fermi LAT on 
the combined analysis of the 10 dSphs by assigning a like- 
lihood of to all samples in our LHC-only scans that have 
an annihilation cross section larger than the 95% upper 
limit in Fig. 2 of Ref. [1^. A more detailed analysis would 
include the full likelihood function in a more refined way, 
but this is not necessary for the purpose of our study. 
We stress that in order to implement the Fermi LAT 
constraints we make the assumption that the neutralino 
found at the LHC contributes 100% to the non-baryonic 
DM in dwarf galaxies, similarly to the consistency check 



approach of Paper I. The resulting two-dimensional pos- 
terior pdf in the planes {mT^,av) and {fl^h'^ , a^_p) are 
shown in the first column of Fig. |3] (filled contours), where 
they are compared with the case where LHC only future 
contraints are used (blue/emtpy contours). Focusing first 
on the top panel, we notice that the LHC alone (empty 
contours) is not going to be able to identify the correct 
mode in the posterior distribution. The secondary mode, 
at large values of av, corresponds to the case where the 
neutralino is Wino-like. This solution can however be 
ruled out once the Fermi LAT information from dSphs 
is included in the likelihood (filled contours). The re- 
sulting identification of the correct cosmological solution 
in terms of the predicted relic density is confirmed by 
the plot in the second panel, showing how Fermi LAT 
limits can eliminate the mode in the distribution corre- 
sponding to subdominant relic density. The last row of 
Fig.|3]shows the ID marginalized pdf for the relic density, 
comparing the LHC only constraints (empty histogram), 
given in Fig. [TJ with what can be otbained by combining 
LHC with Fermi LAT. We see clearly that the Wino-like 
neutralino solution disappears in this latter case. 

The impact of Fermi LAT data can also be observed 
in the reconstruction of the neutralino mass parameters 
of Fig. [21 where Fermi can rule out the region where 
A/2 < Ml. However, Fermi has little impact in the recon- 
struction of the /X parameter, and therefore the Higgsino 
composition of the neutralino is not well determined. As 
a consequence, the reconstruction of the neutralino relic 
density (bottom panel of Fig. [3]) still displays the char- 
acteristic long tail towards small values of il^oh'^. In 
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practice, even the combined data can only constrain the 
relic density within an order of magnitude or so of the 
true value. 



B. Planck constraints from recombination history 

The most robust constraints on the annihilation cross 
section are perhaps those arising from observations of the 
Cosmic Microwave Background (CMB) radiation (47l - l5l| . 
since they do not depend on poorly known quantities 
such as the DM profile at the Galactic center or in dwarf 
galaxies. 

The annihilation of DM particles around rcdshift ^ 
1000 inevitably affects the processes of recombination 
and reionization, modifying the evolution of the free elec- 
tron fraction Xe and the temperature of baryons. The 
evolution of the ionization fraction satisfies the fol- 
lowing equation: 



dxe 1 

^ (1 + z)H{z) 



[R,{z)-h{z)-Ix{z)], (4) 



where Rg is the standard recombination rate |60l . l6ll |. 
Is the ionization rate by standard sources and the Ix 
term represents a "non-standard" source of ionization; 
in our case it takes into account that, during recombina- 
tion, annihilations of DM particles increase the ionization 
rate both by direct ionization from the ground state, and 
by contributing additional Lyman-a photons. Therefore, 
the ionization rate due to DM annihilations can be writ- 
ten as: 



Ix{,z) = Ixi{z) + Ixa{z), 



(5) 



where Ixi is the ionization rate due to ionizing photons, 
and Ixa the ionization rate due to additional Lyman-a 
photons. The rate of energy release dE{z)/dt per unit of 
volume by a relic self-annihilating DM particle at redshift 
z is given by 



dE . 



av 



(6) 



Here, ^Idai is the DM density parameter and Pc the crit- 
ical density of the Universe today; the parameter / indi- 
cates the fraction of energy which is absorbed overall by 
the gas, under the approximation that the energy absorp- 
tion takes place locally. The latter parameter, together 
with the DM mass and annihilation cross section, define 
the parameter we call pann = fcv/m^. 

Both the non-standard ionization rates Ixi and Ixa 
are related to the energy release rate as follows: 



_Cx^dE_ 
nH{z)Ei dt 



Xc 



(1 - C) dE 
nH{z)Ea dt 



(^) (7) 



where Ei is the average ionization energy per baryon, Ea 
is the difference in binding energy between the Is and 2p 
energy levels of a hydrogen atom, nuiz) is the number 



density of hydrogen nuclei and Xi ~ Xa = (1 ^ 2;e)/3 
are the fractions of energy going to ionization and to 
Lyman-a photons, respectively. A fraction of the energy 
released by annihilating DM particles goes into heating 
of baryonic gas, adding an extra term in the standard 
equation for the evolution of the matter temperature (see 
Eqs. (9) and_(10) of Ref. [13). 

In Ref. [47| the angular power spectrum of CMB 
anisotropics was computed, taking into account the pres- 
ence of DM annihilations and an upper limit of 2.0 x 
10~^m^s~-^kg~^ was derived from WMAP5 data at the 
95% confidence level. This limit was already sufficient 
to exclude models that exhibit a large Sommerfeld en- 
hancement for {av), such as those proposed to fit the 
PAMELA and ATIC results. The upper limit is expected 
to become an order of magnitude more stringent, of order 
1.5xl0"'^m'''s"ikg~\ with upcoming Planck data. In or- 
der to study the impact of such data on LHC constraints 
on the pMSSM parameters, we have added a projected 
upper limit on Pann from Planck, shown as the dashed 
red line in the top central panel of Fig. [31 In doing so, 
we adopted a fixed value / = 0.5 for the fraction of en- 
ergy absorbed by the gas. Although the actual value of / 
depends on the particle physics characteristics of the neu- 
tralino, detailed estimations of / (e.g., Ref. ^M,) show 
that this is a reasonable choice. 

As shown in the central panels of Fig. [31 we find that 
future Planck constraints will complement LHC data in a 
very similar way to those arising from the non-detection 
of gamma-rays from Fermi LAT. The combination of 
LHC and Planck will exclude regions of the parameter 
space with high annihilation cross section and low relic 
density, leading again to the indentification of the correct 
cosmological solution for the relic density within an order 
of magnitude (bottom panel). 



C. Constraints from a CTA-like experiment 

We now move on to analyze the implications that an 
actual detection of excess photons would have on the 
parameter space of SUSY DM discovered at the LHC. 
Cherenkov telescopes represent the future (or at least the 
near future) of gamma-ray experiments. Cherenkov tele- 
scopes detect gamma-rays indirectly, through the detec- 
tion of the electromagnetic shower produced by the inter- 
action of primary gamma-rays with the atmosphere. The 
shape of the image created by the shower in the telescope 
camera allows to discriminate photons from hadrons, and 
it also provides information on the incident gamma-ray 
as, e.g. energy and direction (see Ref. for a recent 
review on the Cherenkov telescopes and their role in DM 
searches). We focus here on the upcoming Cherenkov 
Telescope Array (CTA) [g^. Current plans are to build 
a facility on two sites, one on each hemisphere, with tele- 
scopes of three different sizes, for a total of about 100 
instruments. Such large number of telescopes, combined 
with a large field of view, will allow one to monitor more 
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FIG. 3. 2-Dimensional marginal posterior pdf in the plane (rn^, (cv)) (top row) and (Q,^h? , a^^p) (middle row), including 
simulated future LHC data only (blue/empty contours) and adding Fermi LAT upper limit from the combined analysis of 
dSphs (left column, filled regions). The second columns combines LHC data with the upper limit expected from Planck on the 
reionization of the CMB radiation, while the third column combines LHC data with an hypothetical detection of gamma-rays 
from the Draco dSph obtained with CTA. The bottom row shows the ID marginal pdf for the relic density ^.^h? . The inner 
and outer contour encloses 68% and 95% probability regions, respectively. The best fit is shown with the encircled black cross 
while the true value is given by the yellow/red diamond. The black continous line in the bottom row indicates the pdf of Q.^hF' 
for the LHC-only case. 



sources at once, and also to perform surveys of large por- 
tions of the sky [g^, EH. CTA sensitivity is estimated 
to go down by a factor of 10 with respect to current 
Cherenkov telescopes, and the energ y ra nge will cover 
the interval from 10 GeV to 200 TeV [el [11. The CTA 
collaboration will complete its Design Phase in two years 
from now. Currently different prototypes of the final tele- 
scopes are already under construction and will be tested. 



A realistic timescale for the construction of the experi- 
ment is 5 years from now. 

For our purposes, we will consider a rather optimistic 
energy threshold of 20 GeV. Our benchmark point cor- 
responds to a neutralino mass of 139.3 GeV (see Tab. 
so gamma-rays from DM annihilation arc expected to fall 
in an energy range for which a reasonable estimation of 
CTA effective are is 10'*m^ (contrary to the largest values 
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predicted for higher energies). 

We consider CTA observations of Draco, and we as- 
sume for its astrophysical factor an informative prior 
p{X), modeled Gaussian distribution with mean 

1.20 X 10^^ GeV cm^^sr and standard deviation of 
0.31 X lO^^GeV cm-^sr, following Tab. 4 of Ref. (fi^. 
With such prescriptions (and for an astrophysical factor 
equal to the mean of the Gaussian prior), our bench- 
mark SUSY point predicts a flux of A^cta ~ 7.1 photons 
from the halo of Draco above 20 GeV, with 1000 hours 
of observation. As we are interested in the constraining 
power of CTA around the true value of the parameters, 
we neglect Poisson noise in realizing the counts and we 
thus assume an observed number of photons above 20 
GeV iVcTA = 7. There is not enough statistics to build a 
gamma-ray spectrum, so we will only consider the infor- 
mation coming from the detection of all photons above 
the energy threshold. This is implemented in the like- 
lihood as an additional experimental constraint of the 
following form: 

Ccta{^) = j p{XN{x.)\NcTA)p{X)dX (8) 

where A'^(x) is the prediceted number of photons above 
threshold for a pMSSM parameters value x and an astro- 
physical factor equal to unity, A is the unknown true as- 
trophysical factor for Draco (which we marginalize over) , 
p(AA^(x)|A'cta) is a Poisson distribution in iVcTA and 
p(A) represents its prior. 

The results of augmenting LHC data with a future 
CTA-like detection can be seen in the last column of 
Fig. 131 With a detection of only 7 photons from Draco, 
the LHC constraints on the DM self-annihilation cross 
section collapse from 6 orders of magnitude down to 
about 1 order of magnitude. As a consequence, the 
relic abundance can again be constrained around its true 
value. We point out that this dramatic improvements 
hold true even when, like we do here, astrophysical un- 
certainties are marginalized over. 



IV. CONCLUSIONS 

We have discussed the complementarity of indirect and 
accelerator DM searches, and we have shown that cur- 
rent upper limits on the DM self-annihilation cross sec- 
tion arising from Fermi LAT constraints on gamma-ray 
from dwarfs are already sufficiently powerful to break de- 
generacies in the phenomenological parameter space of 
DM that would be present even in high-luminosity LHC 
data. Similar results are obtained if we consider the in- 
formation that Planck will provide on the distortion of 



the CMB, or if we assume a possible future detection of 
gamma-rays from Draco by the CTA telescope. 

In Paper I, we had performed a similar analysis for 
a combination of accelerator and direct detection data, 
and demonstrated how a simple yet powerful scaling 
Ansatz (stating that the local contribution of neutrali- 
nos to Galactic DM was the same as on average in the 
Universe) led to the correct identification of the cosmo- 
logical solution. In this paper, we focused on the power of 
indirect detection techniques, which do not require such 
an Ansatz. Instead, we assumed that the neutralino seen 
in colliders makes up the totality of the DM in the Uni- 
verse (this is what in Paper I was called "the consistency 
check" approach). This corresponds to asking the ques- 
tion: can the observed neutralinos make all the DM in the 
local universe? This approach allows one to discard so- 
lutions corresponding to high annihilation cross sections 
without making any assumption on the expansion rate 
of the Universe at freeze-out, therefore bypassing all the 
difficulties arising from all the effects that might possibly 
modify it (see Ref. [g^ for a recent review). 

The removal of the spurious Wino-like mode which 
would survive with LHC measurements alone and the 
identification of the correct cosmological solution can be 
achieved via various indirect detection channels: current 
Fermi LAT limits on the flux from dwarfs are already 
sufficient to this goal (at least for the benchmark point 
studied here), but Planck constraints will give similar re- 
sults when they become available in 2013, while probing 
very different physics. Furthermore, we have shown that 
a detection of gamma-ray from Draco from CTA can lead 
to similarly stringent conclusions, even when astrophysi- 
cal uncertainties are included in the analysis. 

This work therefore makes the case for a vigorous pro- 
gram of indirect searches, especially in the case where the 
LHC finds evidence for New Physics. We have demon- 
strated that in this case indirect detection methods can 
have a crucial role to play, even in the case where the ex- 
periments only return upper limits to the signal. Together 
with the complementarity between LHC and direct de- 
tection methods presented in Paper I, our results thus 
highlight the fundamental importance of a multi-pronged 
approach to DM identification. 
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